Pulsar Kicks With Modified URCA and Electrons in Landau Levels 
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We derive the energy asymmetry given the proto-neutron star during the time when the neutrino 
sphere is near the surface of the proto-neutron star, using the modified URCA process. The electrons 
produced with the anti-neutrinos are in Landau levels due to the strong magnetic field, and this leads 
to asymmetry in the neutrino momentum, and a pulsar kick. The magnetic field must be strong 
enough for a large fraction of the eletrons to be in the lowest Landau level, however, there is no 
direct dependence of our pulsar velocity on the strength of the magnetic field. Our main prediction 
is that the large pulsar kicks start at about 10 s and last for about 10 s, with the corresponding 
neutrinos correlated in the direction of the magnetic field. We predict a pulsar velocity of 1.03 
xlO~''(r/10^''A')'^ km/s, which reaches 1000 km/s if T ~ 9.96 x 10^° K. 
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I. INTRODUCTION 

The creation of neutron stars, often called pulsars, 
through neutrino cooling starting less than a second 
after the collapse of a massive star has long been 
of interest, with a number of processes contributing 
to neutrino production[l|, The strong interaction 
treatment of these processes was refined by Friman 
and Maxwell using perturbative one pion exchange 
and short-range interactions [1]. 

In recent years it has been observed that many 
pulsars move with much greater velocities than other 
stars in our galaxy. This is called the pulsar kick. 
See Ref.[3] for a review. Pulsars with velocities of 
more than 1500 km/s have been observed. There 
has been a great deal of theoretical effort in at- 
tempts to explain the pulsar kicks. For many years 
a number of investigations of asymmetries in the 
hydrodynamics of core collapse have been carried 
out, and in recent years there have been a num- 
ber of reviews 

1,1, Mi- Although some simula- 
tions have found possible pulsar velocities of 1000 
km/s or more, there is no clear proof that one can 
get the observed large pulsar velocities by the initial 
core collapse. There have been several calculations 
of possible asymmetry in the neutrinos produced in 
strong magnetic fields using the URCA process [9] 
and other processes [l3, [HI during the first few sec- 
onds when the neutrinosphere has a radius of about 
40 km. However, the opacities and short mean free 



paths of neutrinos in the neutrino atmosphere re- 
duces the emission, and these pro cesses cannot ac- 
count for the large pulsar kicks 12, Ts*!. There have 
also been calculations of pulsar kicks resulting from 
oscillation to sterile neutrinos [l^. [Tsl. IT^ ITtI . which 
can escape from the neutrino sphere during this early 
period. The recent MiniBooNE experiment ^18j| with 
previous LSND resiults are not consistent with a sin- 
gle sterile neutrino, but allow models with two or 
more sterile neutrinos. 

It has long been recognized that during the later 
period, when the neutrinosphere radius has been re- 
duced to about 10 km, the radius of the protoneu- 
tron star, that the modified URCA process domi- 
nates the cooling of the protoneutron star[l|,i. It is 
also known that protoneutron stars have very large 
magnetic fields. In the presence of such fields the 
electrons produced in the modified URCA process 
will be in Landau levels [Tgl. [20| 

In the present work we derive the asymmetric neu- 
trino emissivity during the period when the neutri- 
nosphere is just within the protoneutron star sur- 
face, and show that due to the electrons being in 
Landau levels pulsar velocities consistent with ob- 
servations are obtained. In the present work the 
contribution from polarization of the nucleons is not 
included. Our work shows that for about a 10 s 
period starting at about 10 s, when the modified 
URCA process dominates and the radius of the neu- 
trinosphere is a fraction of a km less than the neu- 
tron star radius, the main neutrino emission is asym- 
metric, and during this period the high velocity pul- 
sar kicks are generated. One of us presented a pre- 
liminary version of this work at the CosPA 2006 
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Symposium [2l|. Also, the process of temperature 
equlibrium for the electrons in a strong magnetic 
field is being completed [2^ 

Our paper is organized as follows. In Sect. II we 
give the basic quantities in terms of which we cal- 
culate the emissivity. In Sect. Ill we explain the 
important differences between our formulation with 
the electrons in the lowest Landau level giving asym- 
metric momentum emissivity and the previous ones 
with no kick. As explained in Sect. IIIB, we in- 
corporate the Landau wave function in the lepton 
trace and give the results of the calculation of the 
traces, explaining that with the contributing elec- 
trons all moving in the direction of the magnetic 
field the intregration over the direction of the neu- 
trino momentum is not present. In Sect. HID we 
give our results for the asymmetric emissivity that 
gives the pulsar kick. In Sect. IIIDl we evaluate the 
essential ingredients of proto-neutron star structure 
that is needed for obtaining our final numerical re- 
sult, given in Eq. (26) and Fig. 4. In Sect. IV we 
present our conclusions. In the Appendix we give 
details of the calculation of the matrix elements, the 
nucleon and lepton traces, and the angular integrals. 



II. MODIFIED URCA PROCESS IN A 
STRONG MAGNETIC FIELD: LANDAU 
LEVELS 
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for cooling of protoncutron stars has been treated 
in many publications [U, 0, @]- In a detailed 
calculation^] the one pion exchange (OPE) and a 
short-range interaction were used for the nuclear in- 
teraction, illustrated in Fig 1. 

Diagrams (1,2,3) are for neutral and charged pion 
exchange, diagrams (4,5,6) are the exchange dia- 
grams, and 7 is the short-range diagram. It can 
be shown that for our calculation of the asymmetric 
neutrino emissivity the short-range n-n interaction 
is negligible, so we only consider the OPE diagrams 
here. 

The OPE factor is given by the standard nonrel- 
ativistic form, while the weak interaction used for 
the nucleons is the nonrelativiatic form of the stan- 
dard model. As we shall show, only the weak axial 
interaction, Wa, needs to be considered 



VoPE 

Wa 
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T <7i ■ k- 



7^2 ■ k 



(2) 



= *(9")7m(1- 75)^(9'^) , 



with (T, r the Pauli spin and isospin operators, (1,2) 
refer to the two nucleons at pion vertices, G = 
"^^a ; gA — 1.26, the X the nucleon spinors, and 
the lepton wave functions are ^'(g'^), 5" ((?''), where g*^ 
and are the electron and antineutrino momenta, 
respectively, with the notation q'^ = q^,q'^ = q^. 
Another point that should be mentioned is that the 
inverse modified URCA processes which we neglect 
in the present work. 
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will also have the electrons in Landau levels and 
thereby contribute to the overall pulsar kick, as we 
show below. 



Landau Levels and Product Matrix 
Elements 



FIG. 1: Modified URCA diagrams with OPE and a 
short-range n-n interaction 



The modified URCA process, 

n + n — > 71 + p + e'^ -\- v'^ 



(1) 



The neutrino wave function is a Dirac spinor, 
while the electrons in a very strong magnetic field are 
in Landau levels, which are similar to that of a Dirac 
particle in a plane wave state moving in the z or - 
z direction, where z is the direction of the magnetic 
field, and the motion in the transverse direction is of 
the form of a tightly bound state. For a discussion 
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of Landau wave functions see Refs[19j, |20[. These 
states are labeled by a principal quantum number, 
n, and spin and momentum. 

An essential aspect of the present work, is that 
a sizable fraction of the electrons are in the lowest 
Landau level (n=0). One can estimate the fraction 
of electrons in this n=0 level from the energy gap 
between the n=0 and n=l levels, and T. The pro- 
toneutron stars that receive a large velocity from 
the mechanism that we are considering have a large 
magnetic field strength, about 10^^ to 10^^ G at the 
star surface at a time of about 10 sec, when the 
neutrinosphere is near the surface of the protoneu- 
tron star. From expressions given in Ref[20|, the en- 
ergy gap between the n=0 and n=l levels is about 
6.0 MeV. Our work, shown below, concludes that a 
temperature T with kT about 8.59 MeV will yield 
a pulsar velocity of about 1000 km/sec. From this 
we can use basic thermodynamics to estimate that 
the occupation probability of the n=l state is about 
0.5, and that the overall probability of the electron 
being in the n=0 level is about 40%. 

A crucial point is that the electron in the n=0 level 
has its spin in the -z direction, which we show below 
causes all of the emitted neutrinos to be correlated in 
the z direction, while electrons in the higher n levels 
have both helicities, and give no net pulsar kick. 
Since we assume that all electrons produced via the 
modified URCA process are in the n=0 Landau level, 
this introduces an error of a factor of aqpproximately 
2, but as we shall see this is a trivial correction due 
to the very strong dependence of the resulting pulsar 
velocity on T, and the uncertainty in the magnitude 
of T. 

The lowest Landau level, with n=0 has the form 



s 



i + m 

V''^""'""^"=°(<Z1,9^,0) (3) 



where u is a negative helicity Dirac spinor, = 
qf'-if,, and 7 5ml] ^ = B/{2Bc) ~ B/8 x 
W^Gauss ~ 12.5 to 125 for B^IO^^ to 10^^ 
Gauss. Using E+m ~ E, the wave functions have the 
conventional normalization J = 1. The vector 
weak interaction is not included in our calculation 
as its contribution to the neutrino asymmetric emis- 
sivity is much smaller than the axial vector. The 
pulsar velocity is independent of the magnitude of B 
except for the n=0 occupation probability. 



Although the z-component of the electron momen- 
tum can be in the -\-z or -z direction, we shall show 
that there is no neutrino emission for the electron 
moving in the -z direction, due to the vanishing of 
the lepton trace. 

For the calculation of the neutrino emissivity and 
the pulsar momentum one needs the traces of the 
product matrix elements. The axial product matrix 
element, |Afyip, is obtained by taking the nucleon 
traces over the product of the leptonic traces times 
the square of the weak-strong product: 

|M^P = Tr{nucleon)\Tr{lepton)[l\lj\ 

[Wa{Vope + exchange)]^j\^ , (4) 

with VoPE given in Eq^ and the (i,j) indices are the 
spatial components of the lepton currents, defined 
in Eq.([21). This is treated in detail in the following 
section and in the Appendix. 



We introduce the following notation: 



Pl~P3 , P^P2-Pl 



k + p 



A = 



1 



(5) 



(X)2 G „ 1 



P ml 



R{k) 



[k^f +ml ■ 

The factor R{k) is used in exchange matrix elements. 



III. NEUTRINO ASYMMETRIC 
MOMENTUM EMISSIVITY 

With a strong magnetic field the electrons are in 
Landau levels. That is, the motion of the electron 
transverse to the magnetic field is compressed, and 
the electron is essentially a one-dimensional plane 
wave Dirac particle moving in the direction of the 
magnetic field with energies given as Landau state 
energies. With a very strong magnetic field, such as 
near the surface of the protoneutron star when the 
modified URCA process dominates neutrino emis- 
sion, the electron will fall into the lowest Landau 
level, with the wave function given in Eq.Q. The 
polarization of the nucleons is a very small effect [l|. 
The electrons being in Landau levels gives asymme- 
try to the neutrino emissivity. As we shall now dis- 
cuss, the standard formulation for obtaining the en- 
ergy emissivity from the modified URCA processes 
leads to the momentum of the emitted neutrinos cor- 
related in the B-field direction, and thus directly to 
the pulsar velocity. 
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A. Neutrino Asymmetric Emissivity 

The neutrino emissivity is given in general form 
in many papers, e.g., see Ref|js|: 

^ (27r)3 2tj'^(27r)3 J (27r)3 

(2^)' E ^^'-^A^A^^i (6) 

^{^^ final ^ initial) final Pinitial) 

\Ma? - Mlz^ + M^^ + M^I^ + M^ZI , 

where J- is the product of the initial and final Fermi- 
Dirac functions corresponding to the temperature 
and density of the medium. We use the terminol- 
ogy that Ma is the axial matrix element, which is 
given as the sum over the six Feynman-like diagrams 
shown in Fig. 1, but keeping only the axial weak in- 
teraction, shown in Eq.((2l); stands for the direct 
diagrams 1,2, and 3 in Fig. 1, for diagrams 4,5, 
and 6, M^rj' is the product of (A/j')^ and . 

M^zf and Afj^jJ, M^Za^ are analogously defined. 
The sum runs over the nucleon spins, Si, and the 
neutrino spin. The electron is in the lowest Landau 
level, with its wave function given in EqQ. 

As in Refs the nucleons and the electrons 

are in thermal equilibrium. In this present work 
we assume that the proton quickly reaches thermal 
equilibrium, while the process of the electron state 
transforming to the n=0 Landau state does not in- 
terfere with the proton reaching its Fermi momen- 
tum. Therefore we can use the values for the mag- 
nitudes of the nucleon and lepton momenta derived 
in 

In the present paper we neglect the polarization of 
the nucleons, so the A^ = 1/2 (see Eq.([S])). Thus the 
entire asymmetry of neutrino emission, which causes 
the pulsar kick in the present work, arises from the 
electrons being in the lowest Landau level and the 
modified URCA process. 



B. Neutrino asymmetric momentum emissivity 

As seen from Eq.([S]) our calculation of the asym- 
metric neutrino emissivity starts with the standard 
theory, except that the electron is in the lowest Lan- 
dau level, rather than in the usual Dirac state. This 
turns out to be the crucial point. With the elec- 



tron in the lowest Landau level, the energy emissiv- 
ity is also the projection of the momentum emissivity 
along the B axis. That is, in the following sections of 
the present paper we shall show that the emissivity 
given by Eq.® is of the form 

e''{q^y = angular and energy integrals X 171^1^(9")^ 
= P^^(q'^)^c = -p(ns)c, (7) 

where Tr{lepton)Tr{nudeon)\MA\'^ = \M\'^{q'')Z 
c= speed of light, P^^{q'^Y is the projection of the 
momentum emissivity (momentum/ volume/time) 
along the z=B direction, and p{ns) is the momen- 
tum/volume/time of the proto-neutron star. It 
should be noted that the form equivalent to that of 
Eqll?]) has been used by a number of authors study- 
ing pulsar kicks, such as RefsflGl. IT7|. 

An essential ingredient in our framework is that 
electrons in the lowest Landau level moving in the 
direction of the magnetic field are emitted while 
those moving in the opposite direction are not, 
which produces the asymmetric momentum emissiv- 
ity. This arises from the lepton traces (see Eq.Q 
and Eas. (|38l39p ) It is a straightforwrd exercise in 
Dirac algebra to show that for an electron in the 
lowest Landau level (/ — 'y5)u~ {q^)u~ {q^)] vanishes 
ior q"^ ■ B = —1, so that we can write 

J d^qlTr[lll,] ~ &7:E'[{qy5a + iq''y5,, 

-~S,^{q-^)^]{f = B = z) , 

which is Eq. (PD|) . derived in the Appendix. Note 
{q^y is the ith component of q", and {q")^ = {q^Y ■ 
This shows that when deriving the momentum asym- 
metry we can assume that the electrons in the lowest 
Landau levels are moving in the direction of the mag- 
netic field. We shall take q*^ in the z (B) direction in 
the remainder of this paper. 

An important result of our work is that the pulsar 
velocity does not depend directly on the strength of 
the magnetic field, as is indicated in Ea ipO]) . This 
can be seen from the fact that the integral over the 
transverse components ofcthe electron in the lowest 
landau level, J dg^q^ g-l?!) /2A _ 2A, gives a factor 
of A which cancels the factor of 1/A from the square 
of the electron wave function (see Eq. ([3]) . Therefore 
the only dependence of our pulsar velocity on the 
magnetic field strength is the requirement that it is 
large enough that a sizable fraction of the electrons 
are in the n=0 Landau level. 

Therefore, with the electrons produced with the 
neutrinos in the modified URCA process undergo- 



5 



ing a transition to the lowest Landau level in the 
very strong magnetic field near the protoneutron 
star, the standard theory for energy emissivity also 
gives the emitted momentum of the neutrinos, and 
therefore the recoil velocity of the resulting neutron 
star. However, there is a great deal to discuss, in- 
cluding the calculation of the various traces and the 
effective volume and time for the modified URCA 
process to give a pulsar kick. 

We start with the calculation of the traces of 
|M4p. As an example, let us calculate the direct 
axial matrix element. Writing the product of ma- 
trix elements needed for the emissivity, as shown in 
Eq. dn]) , and using the notation given in Eq. ([5|) , with 
Ml, M2, and M3 the 1, 2, and 3 diagram in Fig. 1, 

M1 + M2 = Axl(T-kxixl[<J-k,(7-l]+X2 
= 2yl? • kxla ■ kxix\x2 
-2AX3CT • fccr • lxix\(^ ■ kx2 (8) 
M l + M2 + M3 = 2AI ■ kxla ■ kxi 



M3 = 

M? ^ 



X4X2 



X3CT • ka ■ Ixixlo- ■ kx2 ■ 



Further results and the calculation of the M 



D-D 
A-A ' 



M 



D-D 
A-A ' 



and M 



D-E 
A-A 



traces of the product matrix 



elements are given in the Appendix. It is important 
to note that the Landau wave function is contained 
in the lepton current, h. 

From the Appendix, Eqs. (l32|37l40|) . and defining 

MSI - + 



M 



E-E 
A-A 

DE 



dWMlzI = j d\lTr[l%]A^{ehk,+k^5,,) 

= SnA^Eek^i-iq^fk^ + 2{q''Yhk^) 
R{kfM2zI{k r) (9) 

- j d?qlTr[l\lj]A^R{k) 

5 2 

— ' ^ i^^i^j ) '2ik,j^kjk 

+2hkj{k^ f - (A: X x 
+(fc • k^f5,.j\ 

~87TA^R{k)Ee[~5k ■ k^ik^k ■ q" 
+k^k'' ■ q") ~ 2q ■ k X k^ik X k")^ 
+Ak ■ q^KXk^ f -f Ak" ■ q^klikf 
+i{q''y-{{k-k''Y ~P{k^f)] . 



C. Angular integrals with k and k*^ angles 
independent 

Following the prescriptions of RefQ the magni- 
tudes of the momenta are given by the Fermi mo- 
menta, so one only does integrals over the ener- 
gies and angles of the momenta, and the momentum 
transfer of the pion in the direct term, k = pi — 
is introduced as an independent vector by inserting 



/ 



d^6ik ^pi+p3) = l. 



For the D-D (direct-direct) term, the only angular 
integral is over the direction of k. For the E-E and 
D-E terms, however, one must deal with k'^ — k + p 
(see Eq.®). Using the (5-functions, one can see that 
neither k'^ nor p are completely independent of fc, 
but with the proton momentum being smaller than 
the neutron, following arguments in Rcf[3], it is a 
good approximation to assume that k'^ and k are 
independent. 

Using the results given in Eq.® for the traces 
of the products of the matrix elements, the integrals 
given in Eq. (|4ip , and the approximations k = k'^ and 
R{k) = 1, with the notation / / = angular integrals 
over k, k'^, and one finds 



M 



D-D 
A-A 

DE 
AA 



Mi 



-E 
-A 



128 
~3~ 



M 



32 



TT^A^k^q'') 



TT^A^k^q-^yq" 



(10) 



Recognizing that the nucleons are in thermal equi- 
librium, the magnitudes of the nucleon momenta are 
given by the Fermi momentum, pp, and therefore the 
integrals over the magnitude of the four nucleon mo- 
menta are carried out via delta functions. Dropping 
the E'^ terms, this gives the final result for the traces 
and integrals over the axial product matrix element 



IM^P = -0.81 X 10^ A^q^iq^ypj, 



-0.81 X 103(X)4^ 



Pp + mi 



niT^ 2llP- 



The next step in the calculation is to carry out 
the nucleon angular integrals. 
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D. Neutrino asymmetric emissivity 

From Eq. ([6]) , based on the general ideas of Ref f3l| , 
but with our formaUsm for the phase space integrals, 
with the angular integrations given by Eqs. (|10|lip in 
subsection III.C, we obtain the neutrino asymmetric 
emissivity: 



\Ma\ 



(27r)9 



PF{e)I . (12) 



The energy integrals, I, are the same as those in 
Refill; and with q" = 85MeV, and q" = A.TkT, we 
find: 



/ = 9.04 X lO^(fcT)^ 
From Eqs. (|llll2ll3|) we obtain 



(13) 



AS 



T 



0.64 X 10^^^^) Vg citC^s 
Pnsc(volume^"'^ time^"'^) , 



3 „-i 



(14) 



where Pns is the momentum given to the neutron 
star, the volume is the active region for the modi- 
fied URCA process and the time represents the time 
interval, pns has the magnitude of e^^ /c in the op- 
posite direction of the net neutrino momentum. 

Due to the short mean free path of neutrinos 
within the neutrinosphere, the main asymmetric 
emission from the process we have proposed will take 
place in the volume between the neutrinosphere and 
the protoneutronstar surface during the period when 
the neutrinosphere is just within the protoneutron- 
star surface. At this time the temperature is ex- 
pected to be in the range lO^K < T < lO^^K for 
a period of 10s starting at about 10s. Taking the 
neutrinosphere and protoneutron star to have radii 
R'^ and Rns, respectively, in km units, we find for 
the momentum given to the pulsar for this period of 
10-20 s 



Pn 



0.43 X 102^( 



(i?'')^) gmcms-\ (15) 



where we use the effective volume for neutrino emis- 
sion as Vns - K-sp/iere = f ^ {RL -(R^f); whcrC 

f is the fraction of neutrinos which escape without 
striking the neutrinosphere for various positions of 
emission. A rough quess is 0.5, and by integrating 
g^, where z is defined in the direction of B, wc find 
that f=0.52. 



1. Radius of neutrinosphere during modified URCA 
emission 

The final step in our derivation is to estimate the 
volume in which the neutrino emission takes place 
with the modified URCA process in a strong magntic 
field. Referring to Ea. ((23)) . we must find the radius 
of the neutrinosphere, i?", assuming that the radius 
of the protoneutron star at this time is 10 km. 

To do this we use the Spherical Eddington model, 
which has been used by a number of authors to study 
the neutrino atmosphere associated with the cre- 
ation of a pulsar [11 |l3,[2l[2l. We follow the recent 
method of Barkovich et. al.[l3l to find the neutri- 
nosphere radius, and estimate the time and temper- 
ature during which our process is taking place. See 
Ref [1^ for a review of the evolution of the birth of a 
neutron star. 

Our starting point is the energy-momentum ten- 
sor, r^'' for the neutrinos with a distribution func- 
tion f'^{x,k,t) for each type of neutrino, giving an 
energy density, U and momentum density F 



U 

pi 



rpoa 



(2^ 
(fk 
(2^ 



kor 
k^r . 



(16) 



Making use of the Boltzmann equation, and recog- 
nizing that the neutrinos have a very short mean free 
path, X'^ , the neutrino distribution can be written in 
terms of the equilibrium distribution, {f^''Y{x, k) as 

r ^ [f'^Y ~Yk- /^if^^Y 

which form we have used to obtain the integrals in 
EqdlSl). 

The mean free path for each neutrino type in a 
medium of nucleons with density p can be written 
as 



1 



X^k'^p 



(18) 



with the constant determined from the constants 
of the standard weak interaction model and the cross 
sections for each neutrino type. Assuming the spher- 
ical Eddington model one finds for the energy and 
momentum densities, including the electron, muon, 
and tau neutrinos, for equilibrium temperature T 



7 



U 
F 



111 

40 



rpl 

1 dT^ 



36xp dr 



(19) 



The time dependence of T can be found from 
energy-momentum conservation: 



Quj^i^v ^ oor 
af/ + A-i? = 0. 



(20) 



Using the equations of hydrostatic equihbrium, 
Barkovich et. al. showed that the T can be obtained 
from the first order differential equation 



rfT - a{x)) 

d^^ 'T22;2(l -a(a;)) 



= 



(21) 



with T = T/Tc, X — r/Rns and a{x) depends on 
the constant be and two other constants, as well as 
T. We use the constants of RefflTl, with = 40 
MeV, slightly altered for our expected luminoscity in 
the 10-30s interval when the modified URCA process 
dominates. 

Our solutions are shown in figures 2 and 3. 
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FIG. 2: 
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FIG. 3: 



In our solutions we find that a(x) << l,a{x) << 

_ - 2 

T,a{x) ~ T{x) {x > 1), and that the result for the 



mean free path is 

(A^)-i 



T(T2 - a)cTO-i 



(22) 



With this solution the radius of the neutrinosphere 
can be obtained in terms of R„^ from the relation 



R''/Rr^ 



T{T^ - a)dx 



2 cm 

3 Rns 



(23) 



Using our solutions we find that for Rns =10 km, 

R" ~ 9.96 km (24) 

when the temperature is in the range T ~ IG^^if , so 
that from Eq ()23p the neutron star momentumm is 



5.14 X lO^'^gm cm/s(Y^)'^ 



(25) 



For a neutron star with the mass of the sun = 2 x 10'^^ 
gm, including a factor of 0.4 for the n=0 occupation 
probability, 



= 1.03 X 10-'*(-— )^km s-i 



101 



(26) 



giving a velocity of v~ 1000 km/s for T~ 9.96 x 10^° 
K, which is in the expected range. Fig 4 illustrates 
the velocity of the pulsar as a function of T. 

Therefore we find that the modified URCA pro- 
cess can produce the observed velocities of 1000 
km/s if T or more during this period if the temper- 
ature is sufficiently high when the neutrinosphere 
is slightly within the protoneutron star due to the 
electrons being in Landau levels. These large pulsar 
kicks start about 10 s after the supernova collapse. 



6.0 S.O 10.0 12.0 14.0 

T ( lOlO K) 

Pulsar velocity v.s T assuming constant emission volume 

FIG. 4: 
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IV. CONCLUSIONS 

The basis of our present work is that for a mag- 
netic field strength strong enough for a sizable fac- 
tion of electrons produced in the modified URCA 
process to be in the lowest Landau lavel, all of the 
neutrinos emitted will be correlated with the direc- 
tion of the magnetic field, giving a pulsar kick. The 
resulting velocity depends strongly on the tempera- 
ture, and only indirectly on the magnetic field. Us- 
ing expected values of 10^^ to 10^^ Gauss for B, and 
the T for which we get a pulsar velocity of 1000 
km/s, we estimate that about 40% of the electrons 
are in the n=0 Landau level, which is satisfactory. 

The calculation of the asymmetric neutrino emis- 
sivity, and the resulting pulsar velocity arising from 
electrons being in Landau orbits in the strong mag- 
netic field near the protoneutron star, is straight- 
forward with the modified URCA process. Using 
the standard properties of protoneutron stars in the 
10-20 s time interval, we find that with the elec- 
trons in Landau levels, the modified URCA pro- 
cess can account for the measured pulsar kicks for 
T > 9.96 X 10^"K, with the neutrinosphere surface 
just inside the protoneutron star surface. Studies 
give T ~ lO^^K, or even greater near the protoneu- 
tron star surface. We predict a strong correlation 
between the protoneutronstar T and pulsar velocity, 
as well as a strong correlation between the direction 
of the pulsar's velocity and the direction of B. Since 
the luminoscity of the pulsar is related to proper- 
ties of the protoneutron star, such as strengths of B 
and T, one can also expect our result to predict a 
correlation between v and L of the pulsar for high 
L. 




In figure 5 the correlation between large L and 
large v that has been observed is shown. It is a 
subject for future study for us to see if our model is 
consistent with that observation. 

Our unique prediction is that the main neutrino 
emission during this period when the neutrinosphere 
is just inside the the neutron star is almost entirely 
asymmetric, when the pulsar kicks should occur, 
with a strong correlation with Ei,. E.g., E^, = 35 
MeV^ kT = 11 MeV and v==2,000 km/s. 

It is interesting to consider the neutrinos ob- 
served in the Kamiokanda-Il[2^ and IBM^3] detec- 
tors from supernova SN1987A. In the analysis of the 
neutrino data[2^, in which the Kamiokanda-II and 
IBM data are plotted as a function of time, there 
might be a gap in the data, with additional neu- 
trinos seen after 10 s (with the first neutrinos seen 
at about 1 s). See Fig. 6 for KAM II and IBM 
results }2a]. 

Although our agreement with the time gap seen 
in Fig 6 is not statistically significant, we predict 
asymmetric neutrinos start to appear at about 10 
s, correlated in the B direction. This should be ob- 
servable in future measurements of neutrinos from 
supernovae with today's improved detectors, with 
neutrino energies determining the pulsar velocities. 
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FIG. 6: Energies of neutrinos from SN1987A 



FIG. 5: Pulsar velocity vs luminoscity, RefQ] 
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V. APPENDIX 

A. Axial-Axial (A-A) Matrix elements 

The axial OPE matrix element has a direct and 
an exchange part, 



Ma 



MS 



Ml, 



(27) 



corresponding the diagrams 1,2, and 3 for the direct 
and 4,5, and 6 for the exchange in Fig. 1. Using the 
notation Xn for the spinor of the nucleon n (n=l,2 
for the initial neutrons, n=3 for the final neutron, 
and n=4 for the produced proton), from Ea. (|2l6p 



M 



A = 2A[l ■ kxl<7 ■ kxixlx2 - xlo- ■ ka ■ Ixi 



xlo- ■ kx2 



(28) 



IS 



Mf = ^R{k)MS{k ^ r , 1 ^ 2) (29) 
Using Xix\ = the D-D product matrix element 

Ma- A = ^A^Tr[l ■ kl'^ ■ k{Al)[A2) + {A?,){A'i) 
-l^ ■ k{A5)(A6) - I ■ k{A7){A8)] , (30) 



where, assuming that k = k — q 
with the notation k k 



k, A„ = 1/2 and 



(Al) 
iA2) 

m 

{A6) 
(AS) 



Tr[a ■ fcAgCT • Mi] = -fc^ 
Tr[cr ■ l^a ■ kA^a ■ ka ■ lAi] 



-Lljk 6ij 



2" 



{lA) 

Tr[(T ■ fcA2A4] = 
Tr[(j ■ fcA4A2] = 



(31) 



From this we find 



M^Ia = A^Tr[lllj]k\hk, + kH,,) 



(32) 



MfzS = A'R{krTr[lll,]{knktk^ + ikrS^,) 

Using k ~ kjk'^ ~ fc'^, the two direct-exchange 
product matrix elements can be written as (note 



A-A 



M'izXik^k^)) 



M 



D-E 
A-A 



= -4A^R{k) V [l^ ■ kl ■ k^DEl (33) 



spins 



+DE2 - ?t . kDE3 - I ■ k^DEA] , 



with (assuming A„ ~ I /2 for all nuclcons) 
DEI 



—Trla-ka-k"] = ifc • F 
16 ^ ^8 

DE2 = — Trier ■ l'^ a ■ ka ■ k^a ■ la ■ ka ■ k" 

16 ^ 



"o^l^it ^ ' ^ {kii j -|- k^ij^ -\- kikj(^k ) 



e\2-\ 



+kl^{kY - (fc X r),(fc X r)j + %(fc X k") 



DE3 



DEA 



-Tr\a ■ ka ■ k^a ■ la ■ k" 



16 

-{2k ■ kn ■ k" - I ■ kik^f) 
8 

— Tr\a ■ l^a ■ ka ■ k'^a ■ k] 
16 ^ ^ 



= -{2k-kn'' ~l^ ■k'-'ikf) . (34) 
8 

To obtain Ea. ([M)) we use the trace relationships 

Tr[a ■ Aa ■ Ba ■ Ca ■ D] = 2{A-BC ■ D (35) 

-Ax B -C X D) , 



Tr[a ■ Aa ■ Ba ■ Ca ■ Da ■ Ea ■ F] 
= 2E ■ F{A ■ BC ■ D - A X B ■ C X D) (36) 
-2[A -BExF-CxD + C-DExF-AxB 
~{{A X B)x {C X D)) ■ E X F] . 

From this, defining Mff = M^Ia + M^Za we 
obtain 

M|f = -A^R{k)Tr[l\l,][-h-k^{h^ + klk,) 
+2hk,{k' f + 2km{kf + (fc • k^ f6^, 



-{k X k'')i{k X fe^)j 



B. Lepton Traces 



(37) 



The lepton trace, Tr[4^j], from Eq.®, with the 
neutrino in a standard Dirac state, u''^{q'^), and 
the electron in the lowest Landau level, 'i'{q'^) = 
i(^/7)~-'^)e~'*^^ /'^^'''^w~(q'^), with the electron mo- 
mentum along direction of the magnetic field, its 
spin in the opposite direction, is defined as 

Triljl,] = le-^'il)y^Tr[j.fj^{I-j,) 

u-{q'^)u-{q'^)] . (38) 
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Using me ^ E^, and dropping terms in {q^Y, which 
do not give rise to momentum asymmetry, 



7» fljU^j 
{I-j5)u-{q')u~{q^) 



f - '^{qrir (39) 



[7° + 7°7^ + 7^7^ 



W2 1 

a q± — e 
7 



{r = B^ z) 

M.?/! ^ 27r. 



An essential part of our work is that (/ — 
{q'^)u^ {q*^) vanishes if the electron momentum 
is opposite to the B-field, and that only electrons in 
the lowest Landau level in the direction of the B- 
field contribute to the emissivity. From Eqs. ()38|39p 
we find for the electron in the lowest Landau level 



d?qlTr[l\lj\ ~ S^E^q'^yda + iq^r^i^ 



5.,{qn'W ^ B ^ z) .(40) 



C. Angular Integration 

Using k = pi — p3,fc^ = P2 — Pi as indepen- 
dent variables (see text for discussion), and defining 
J J = J dilk J We need the following angular 
integrals 



ki — 



K = (41) 



k-Ak-B = ^A-B 



k- Ak- Bk-Ck- D 
A-CB-D + A-DB-C) 



{k-k^f 



i{kxk'),f 
k ' k k ^k ' (J 



^^iA-BC-D- 
15 



47rfc^ 

^k^ 
3 



47r 



k\k 



^k'ik^r 

k\k^fq. 



47r 



klk'' ■ q 
fc X fc"^ • q{k X k'^')z 



-^k\k^)\. 



(k X fc'=),fc -k^ = . 
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